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Abstract 

In two models it is shown that a hght pulse propagates from a vacuum 
into certain media with velocity greater than that of a light in a vacuum (c). 
By numerical calculation the propagating properties of such a light are given. 



Recently L. J. Wang and his collaborators in their experiment [1] have found 
that the group velocity of a laser pulse in particularly prepared atomic caesium gas 
can much exceed that of light in a vacuum (c). Here we will explore the problem 
from a theoretical piont of view. Many years ago A. Sommerfeld rigorously proved 
[2] that the velocity of light pulse can not exceed c in absorption media. We call this 
conclusion as basic theorem afterwardz. In the following we will first examine why 
usually the basic theorem holds and in what condition it will be violated. Then two 
models are proposed, where with precise numerical calculation we will show that 
the basic theorem is indeed violated. Therefore the properties of light propagating 
in media of the models are given. We believe that these properties may appears in 
more realistic models. 
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Let a light pulse propagating along x axis toward its positive direction in a 
vacuum (a; < 0), at time t = arriving at x = 0, and then entering into a medium 
{x > 0) afterwards. At x = 0, the amplitude of the light pulse changes with time t 
as 

/w4^"^ (1) 

^ ^ \ t<0 . ^ ^ 

Usually f(t) can be rigorously expressed in Fourier integration as 

f(t) = Re / A(n)e~''''dn . (2) 
Jo 

For simplicity suppose there is no reflection of light at x = 0. The amplitude of the 
light pulse entering into the medium (a; > 0) is [2] 

roo 

f(t,x)^Re A(n)e-^"*+^*^^dn , (3) 
Jo 

where k = nfi{n)/c and fi{n) is the complex refractive index of the medium, which 
depends on the frequency of incident light (dispersion). In vacuum [x < 0), the 
amplitude of the light pulse is 

g{t,x) = i?ey~A(n)e-^"*+^"^/^dn {x < 0) 

= /(*-x/c)={^('7/'' ;-f ^; (4, 

[ t - x/c < . 

The shape of the light pulse propagating in a vacuum does not change because all 
its Fourier components in (4) propagate with the same velocity c and do not decay. 
In particular, all these components completely cancel each other in the space-time 
region t — x/c<0{oic9 — tc/x< 1) as long as the hght pulse propagates in a 
vacuum. One may think that such a cancellation may not occur when light pulse 
propagates in a medium because of dispersion and absorption. But Sommerfeld 
proved that such cancellation also occurs when light pulse propagates in absorption 
media, which is just the basic theorem. In his proving Sommerfeld let 

, , I sin u t t > Q 
[ t <0 , 

so that (3) becomes [2] 

fit, x) = ^Re I e-^("*-"-(")}/l,) dn , (6) 
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where ii{n) is taken as 



a 
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which is the Lorentz-Lorenz refraction formula. a^,no,p are the constants of the 
medium, no, p represent the characteristic absorption frequency, damping constant 
of the medium. Usually p > 0, light propagating in the medium decays. The 
integration path c in (6) is shown in Fig.l, which is along the real axis of n from 
+00 to — oo through n — u by a small semicircle in the upper half of the complex 
plane. fj,{n) in (7) has branch pionts: 



Ui^2 — —ip^ynQ — p^, where p, — oo 



Ni^2 — —ip i \J -\- — p"^ , where p — . 

We joing Ui to A'^i and U2 to N2 by two branch lines, which lie in the lower half of 
the complex plane, when p > 0. Because there is no singularity and branch lines of 
the integrand of (6) in the upper half plane, one can replace the integration path 
c by u in Fig.l, which is parallel to the real axis in the upper half plane. When u 
moves to infinity in the upper half plane, /i — > 1 and when t — x/c < 0, then (6) 
= 0, which is just the basic theorem. But when p < (this is our model 1), the 
branch lines UiNi and U2N2 lie in the upper half plane and the integration path u 
is not equivalent to the path c in (6). Now the equivalent path u should be taken as 
U1 + U2 + U3 + U4 in Fig.l. By the same argument above, integration along ui vanishs. 
Integrations along a pair of U2 cancel each other. The remaining integations along 
U3 and U4 usually do not vanish because the branch lines lie in them. Therefore 
when p < (model 1), the basic theorem does not hold. In the following we will do 
numerical integration of (6) to show that the basic theorem is indeed violated and to 
see what happens, p < means propagation of light in the medium is gain-assisted 
light propagation. 
From (6) we get 

f{t,x) = - i?e[ze^'^(^)] + —Re / [e^«'(--+^) - eT«'(-+^)]— , (8) 
2 2n Jo z 

where 7^(71) = —int+inpx/c; 7 = xtiq/c, n — n/riQ, D — v/uq, — a? /uq, p — p/rio, 

all of them are dimensionless. Let 

w{z)^X{z)+iY{z) . (9) 
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When z » 1, 

X{z) ~ -a^p/z'^ , Y{z) - z{l -9)- a^/{2z) , (10) 

where 9 — ct/x. 

In reference [2], the typical values of parameters are given as 

no = 4 X 10^^ s'S = 1.24 n^, p = 0.07 no , (11) 

where the medium is solid or liquid. For gas is about 1.001. 

The numerical integration for (8) is difficult when | jY(z) \ becomes very large 
and hence e'''^^^^ is a very fast oscillatory function of z. In fact when no — Ax 10^^ 
and X = 1cm, J = | x 10^ is very large. But if 7 < 100 for example, i.e., x < 
7.5 X lO^^cm, we can do numerical integration of (8) with high precision. Now the 
integration of (8) is divided into sum of /(f and J^. With R ~ 100 for example, 
fast oscillatory integrand appears in J^, when ^ >> 1. Using (10) 

Re{e^^^'^) = e-^"'^/^' \cos[^z{l - 9)]cos^ + sin[-fz{l - 9)]stn^\ , 

I ZiZ ZiZ I 

where the fast oscillatory factors are seperated in forms of sin and cos functions. 
Mathematica can do that kind of integration with high precision. As an example, 
let = 1.24, p = 0.07, P = 10,7 = 1,^ = 0.98 in (8), one may get f{t,x) = 
—1.02057 X 10~^^, which should vanish exactly due to the basic theorem (p > 0, ^ < 
1). Considering the amplidude of incident hght pulse is 1, 10~^^ is a very high 
accuracy of calculation. Now let p — —0.07 and other parameters unchanged, fit, x) 
in (8) is 0.0630255, which is a definite evidence to show that the basic theorem is 
violated indeed when p < 0. Because 7 = xuqIc, 9 = tc/x, fit, x) may be looked as 
a function of 7 and 9, hij, 9). Let us fix 7 = 1, i.e., x = 0.75 x 10~^cm, and 9 change 
from (-7) to 4, i.e., the time t from —7/no to 4/no, the amplitude of light pulse /i(l, 9) 
in the medium is shown in Fig.2, where = 1.24, p = —0.07,7 = l,i^ = 10(z/ = 
10 no). Again when p < 0, /i(l, 9) does not vanish and the basic theorem is violated. 
In Fig.2 even if t < 0, the amplitude hil,9) still does not vanish, which means that 
before the incident light pulse arrives at the medium, the light in the medium is 
already produced. For convenience, we call the light produced in the medium when 
^ < 1 as fastlight and that when ^ > 1 as normal light. When the basic theorem 
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holds, the fasthght vansishes. Some maximal values of light amplitude near 6 — 1 
and their corresponding 9 values are listed in table 1, where two characteristics are 
shown: 

(1) The amphtude of fastlight is oscillatory and decays as ^ — > — oo; 

(2) The period for each oscillation for 9 less and near 1 are roughly around 
^ = 6, which means its frequency is near the characteristic frequency Uq of the 
medium (the corresponding period is ^ = 27r). This near equality is due to the fact 
that the Fourier components of the light pulse having freqencies equal to or near no 
are most gain-assisted when p < 0. When 9 > 1, the normal light soon oscillates 
with the frequency u of the incident light pulse (corresponding period 9 — 0.27r) and 
its amplitude is a little bit lager than 1(1.00089). When 9 — 1 with other parameters 

— 1.24:, p — —0.07,7 — 1 unchanged the two properties (1) and (2) remain, but 
the amplitude of the fasthght increases to 6.14 near 9 — 1. The amphtude of normal 
light inceases to 7.31 with frequency u — uq. 

Model 2: the p in (7) depends on frequency n as 



which is not continuous and therefore p(n) and the integrand in (6) are not the 
analystic functions of n. We still can use numerical calculation to see whether the 
basic theorem holds: 

(a) = 1.24, pi = -0.07, p2 = 0.07,^ = 1,6 = 0.01, P = 1,7 = 1. Now light 
amplitude in the medium are gain-assisted when (1 — 6) no < n < (1-1-6) no and decays 
otherwise. Fastlight with above two properties appears again and its amplitude near 
^ = 1 is 3.81, while that of normal hght is about 4.5; 

(b) pi = 0.02 with other parameters unchanged as in (a). Now light in the 
whole frequency range in the medium decays. Still fastlight remains with the two 
properties, but its amplitude becomes small (0.0757) near 9=1. 

So we may conclude that if p(n) in (6) is an analystic function of n with sin- 
gularities (such as poles or branch lines) appearing in the upper half plane of n, or 
p(n) is not an analystic function of n at all, the basic theorem in general may not 
hold and fastlight appears. The Fourier components of a light pulse now will not 
cancel each other in the medium in the space-time region t — x/c < 0. This is why 
fastlight appears. 




(12) 
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How to measure the velocity of the Ught pulse when fastlight appears in the 
medium? Suppose a light pulse produced in a source at t = 0, propagating a 
distance ii in a vacuum, then going through a medium with thickness £2- Just 
after the medium a light pulse detector is put. Usually the fastlight appears in 
the medium after the light pulse is produced and the amplitude of the fastlight is 
inceasing when the light pulse is approaching to the medium. If the amplitude of 
fastlight is able to become large enough to trigger the detector at t — ti, one may 
take V — •^2)/^! as the velocity of the hght pulse propagating from the source 
to the detector, which certainly exceeds c. 

Although the models proposed above are not completely realistic, we believe that 
production of the fastlight and some its properties in the models may remain in a 
more realistic model, which is under investigation now. 

We would like to thank professors Gu Yi-fan and Dong Fang-xiao for their ben- 
eficial discussion. 
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Table 1. Some maximal values of light amplitude h{l,9) near 9 = 1 
(a^ = 1.24, p = -0.07,7 = 1 and = 10 in model 1) 
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-46.86 


-40.70 


-34.57 


-28.44 


-22.35 


-16.29 


-10.27 


h{i,e) 


0.00909 


0.0125 


0.0172 


0.0232 


0.0316 


0.0413 


0.0580 


















-4.59 


0.8714 
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3.03575 
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1.0089 


1.0089 


1.0089 


1.0089 
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Figure 1: The integration paths of (6). 




Figure 2: Figure 2. ^) (a^ = 1.24, p = -0.07, 7 = 1 and P = 10 in model 1.) 
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